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1. Introduction 

At the LHC the observation of hard processes with many jets, gauge bosons and perhaps heavy 
new particles and their subsequent theoretical interpretation is of primary importance. The quanti- 
tative analysis of these processes requires the knowledge of the corresponding cross-sections and 
correlations at next-to-leading order (NLO) accuracy in perturbative QCD. The standard method 
of calculation based on Feynman diagrams becomes very cumbersome for multi-leg processes. In 
gauge theories the conventional Feynman diagram method produces intermediate results which are 
much more complicated then the final answer. The number of Feynman graphs grows very fast 
with the number of external legs. For a tree-level Af-gluon scattering the number of individual 
Feynman graphs is approximately N^ N ^ (within 5% accuracy up to 16 gluons) [jjj]. This decom- 
position generates a large number of terms. With a growing number of external particles it becomes 
a forbidding task to simplify the resultant expression analytically. Recently more numerical tech- 
niques have been developed (see e.g. ref. [Q]). However, because of the stronger than exponential 
growth in the number of Feynman diagrams these brute-force methods become computationally 
very intensive for amplitudes with six or more legs. 

Unitarity methods have been suggested as an alternative, more efficient procedure for loop 
calculations a long time ago [^, @,|5j]. Their use in the context of gauge theories is especially ben- 
eficial [||, 0]. The computing time is governed by the efficiency in computing tree amplitudes and 
by the number of cuts. In these applications, the unitarity cut is four dimensional. This allows 
the use of helicity method and surprisingly simple analytic answers [^] have been derived. The 
four-dimensional unitarity method reconstructs only the so-called cut-constructible part of the am- 
plitude. The remaining rational part is obtained using known properties of the collinear limit . In 
supersymmetric theories, which have improved ultra-violet behaviour, the rational part vanishes. 



Recently, new ideas (for a review see [Bp) on twistors [10], multi-pole cuts (generalized uni 



tarity) [|l~ij], recursion relations [|12|, [13|, 14, 15], unitarity in D-dimension [16, 17] and the use 



of algebraic parametric integration technique [ |18| ] have made the unitarity cut method even more 



promising. It appears that ultimately one can find an efficient algorithm which can be used to 
calculate the full one-loop amplitudes in terms of tree-level amplitudes. 



Here we briefly describe a semi-numerical four dimensional unitarity method Q19[ ] which ex- 
pands the algebraic method of ref. [ [l8| ] by developing a numerical scheme. The numerical algo- 
rithm evaluates only the cut-constructible part of the one-loop amplitudes. 1 



2. Structure of the one loop amplitude 

The generic D-dimensional Af -particle one-loop amplitude (fig. [I]) is given by 2 

m[Pl,P2,---,PN) =J [dl] did 2 ---d N ' 

where p\ represent the momenta flowing into the amplitude, and [dl] = d D l. The numerator struc- 
ture ^V(pi,p2, ■ ■ ■ ,Pn'J) is generated by the particle content and is a function of the inflow mo- 

'This scheme has been extended recently into a semi-numerical £>-dimensional unitarity method |po|]. For other 
promising methods for calculating the rational parts see Refs. jpl[^]. 

2 We restrict our discussion to (color) ordered external legs. The extension to more general cases is straightforward. 
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Figure 1: The generic Appoint loop amplitude. 



menta and the loop momentum. Since the whole amplitude has been put on a common denominator, 
the numerator can also include some propagator factors. The dependence of the amplitude on other 
quantum numbers has been suppressed. The denominator is a product of inverse propagators 



d i = d i (l) = (l + q i f-mf= ( l-qo+Ypi ] -mf , (2.2) 




where the 4-vector qo parameterizes the arbitrariness in the choice of loop momentum. The one- 
loop amplitude in D = 4 — 2s can be decomposed in a basis set of scalar master integrals giving 

£*n(pi,P2,---,Pn) = Y, a h(.PuP2,---,PN)Ih ,+ Kh(PUP%,---,PN)hh 

\<h<N l<h<h<N 

(Pl,P2,---,PN)Ih 

(2.3) 

1 < i i <h < (3 <JV 1 < i i < h < h < k < N 

where the master integrals are given by h v -i M = f[dl] d . }. d . • Analytic expressions for D- 
dimensional master integrals with massless internal lines are reported in ref. [Q]. The correspond- 
ing results for divergent integrals with some massive internal lines are reported in ref. [p3fl. The 
maximum number of master integrals is determined by the dimensionality, D, of space-time; for 
the physical case this gives up to 4-point master integrals. The unitarity cut method is based on 
the study of the analytic structure of the one-loop amplitude. The coefficients are rational func- 
tions of the kinematic variables and will, in general, depend on the dimensional regulator variable 
£ = (4 — D)/2. When all the coefficients of the master integrals are calculated in 4 dimensions 
we obtain the "cut-constructible" part of the amplitude. The remaining "rational part" is generated 
by the omitted tff(e) part of the master integral coefficients. For a numerical procedure we need to 
recast the study of the analytic properties of the unitarity cut amplitudes into an algebraic algorithm 
which can be implemented numerically. In ref. [ 18 ] it was proposed that one focus on the integrand 
of the one-loop amplitude, 

v / , A jV{pi,p2,...,PN\l) , n A , 

■s/n{Pi,P2,---,Pn\1) = did 2 ---d N ' 
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This is a rational function of the loop momentum. We can re-express the rational function in an 
expansion over 4-, 3-, 2- and 1-propagator pole terms. The residues of these pole terms contain 
the master integral coefficients as well as structures (so-called spurious terms) which reside in 
the "trivial" space, the subspace orthogonal to the "physical" space. The physical space is the 
subspace spanned by the external momenta of the corresponding master integral. The spurious 
terms are important as subtraction terms in the determination of lower multiplicity poles. After 



integration over the loop momenta, Eq. ( [2.3D is recovered. This approach transforms the analytic 
unitarity method into the algebraic problem of partial fractioning a multi-pole rational function and 
allows for a numerical implementation. 

3. Parameterization of the loop momentum on the unitarity cuts 

We restrict ourselves to a 4-dimensional space. Given the master integral decomposition of 



Eq. (2.3) we can partial fraction the integrand of any 4-dimensional Af -particle amplitude as 



&N\})= l_ —j — •+ ]_ -j—j \- ]_ , , , + ]_ , , , , • (3.1) 

l<h<N a 'i l<h<h<N a h a h \<h<i 2 <h<N a h a i2 a k \<ii<i 2 <h<k<N a 'i a h a h "U 

To calculate the numerator factors, one evaluates the residues by taking the inverse propagators 
equal to zero. The residue has to be taken by constructing the loop momentum kj...k such that 
di(kj-k) = dj{lij...k) = ■■■ = dk{hj-k) = 0. Then the residue of a function F(l) is given by 



Re Sij ... k [F(l)}= [d i (l)d j (l)---d k (l)F{l) 
The quadruple and triple pole residues are now given by as 



(3.2) 



ij-k 



d m (l) = Res yJU , c m {l) = Resy* L/ N (l) - £ jf^J . ( 3 - 3 ) 

with similar expressions for the double fe !l( - 2 (Z) and single a,-, (Z) pole residues. As an illustration we 
briefly outline how to construct the residue functions for quadruple cuts. 

4. The quadrupole residue 

To calculate the box coefficients we choose the loop momentum l^i such that four inverse 
propagators are equal to zero, 



dijkiilijki) = Resijki \^n{1)J diQiju) = djikju) = dk(hju) = di(l ijk i) = . (4.1) 
We will drop the subscripts on the loop momentum in the following. Because we have to solve the 



unitarity constraints explicitly, we have to choose a specific parameterization. In ref. [19] the van 



Neerven-Vermaseren basis [ J24f j is used which gives a very natural parameterization of the "trivial" 
and the "physical" space in terms of dual vectors constructed from the inflow momenta for a given 
cut type. We can decompose the loop momentum as 

^=V z f + a 1 ^. (4.2) 
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The variable ct\ will be determined such that the unitarity conditions dj = dj = dk = d/ = are 
fulfilled. V4 is a well defined vector in the "physical" space constructed from the three independent 
inflow momenta; n\ is the unit vector of the one-dimensional trivial space. One finds two complex 
solutions 



l^=V^±i^Vi-mfxn 1 } , (4.3) 

which are easily numerically implemented. We note that the four propagators are on-shell and the 
amplitude factorizes for a given intermediate state into 4 tree-level amplitudes ./^°' . The residue 



of the amplitude in Eq. ( |4.1| ) is given in terms of tree amplitudes as 

Res^(«^(Z ± )) = ^ Hlf-p i+u ...,pj--lf)x^ \lf-p j+1 ,..., Pk --l^) 

(l±;p k+l ,...,pr,-l±)x^ \lf;p l+l ,...,pf,-lf) , (4.4) 



where the loop momenta l„ are complex on-shell momenta and there is an implicit sum over all 
states of the cut lines (such as e.g. particle type, color, helicity). The tree-level 3-gluon amplitudes, 



\ are non-zero because the two cut gluons have complex momenta [jllj]. Any remaining 
dependence of the residue d^i on the loop momentum enters through its component in the trivial 
space, djjkiQ.) = d\jki{n\ ■ 1) • The number of powers of the loop momentum / in the numerator 
structure is called the rank of the integral. After integration we find that («i • I) 2 ~ n\ = 1. Thus 
rank one is the maximum rank of a spurious term (which by definition vanishes upon integration 
over I). Hence the most general form of the residue is diju (I) = dijti + diju l-n\ . Using the two 



solutions of the unitarity constraint, Eq. (1.3), we now can determine the two coefficients of the 
residue 

Res, 7/W (Mr(/ + )) + Res yw (<e^(r)) _ Res yw - Res,-,*/ (V w (r)) 



Hjkl — ^ , aijki 



2iJV 2 -mj 



(4.5) 

After the subtracting the quadruple cut contributions from the amplitude we can repeat the proce- 
dure for the triple, double and single cuts. 

5. Numerical results 



As an application in ref. [19] the 4-, 5- and 6-gluon scattering amplitudes at one-loop have 
been recalculated with the new method. The cut-constructible parts of the ordered amplitudes are 
also known analytically making a direct comparison possible. Also, the 6-gluon amplitude was 
numerically evaluated using the integration-by -parts method [||]. To compare with the analytic 
results 100,000 flat phase space events has been generated for the 2 — > [n — 2) gluon scattering. 
The events are required to have cuts in order to avoid soft and collinear regions in the momenta of 
the outgoing gluons. The evaluation time for 10,000 events is: for a 2 — > 2 gluon ordered helicity 
amplitude 9 seconds, for a 2 — > 3 gluon ordered helicity amplitude 35 seconds and for a for a 2 — > 4 
gluon ordered helicity amplitude 107 seconds. Note that using the integration-by-parts method of 
ref. [Q] the evaluation time for 10,000 events would be approximately 90,000 second. The six- 
gluon evaluation is only three times slower than the five gluon evaluation and eleven times slower 
than the four gluon amplitude. 
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6. Summary and Outlook 

The numerical unitarity method provides an efficient method to evaluate next-to-leading order 
corrections to multi-leg hard scattering amplitudes. It is applicable for processes including massive 
and massless particles as well as bosons and fermions. Very recently it has been generalized to 



D-dimensions so one can reconstruct the full amplitude including the rational part [gOp. We expect 
that in the future it will be used in a number of important physics applications. 
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